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ABSTRACT 

Transfer ribonucleic acid (tRNA) modifications, es- 
pecially at the wobble position, are crucial for proper 
and efficient protein translation. MnmE and MnmG 
form a protein complex that is implicated in the 
carboxymethylaminomethyl modification of wobble 
uridine (cmnm^U34) of certain tRNAs. MnmE is a 
G protein activated by dimerization (GAD), and ac- 
tive guanosine-5'-triphosphate (GTP) hydrolysis is 
required forthe tRNA modification to occur. Although 
crystal structures of MnmE and MnmG are available, 
the structure of the MnmE/MnmG complex (MnmEG) 
and the nature of the nucleotide-induced conforma- 
tional changes and their relevance for the tRNA mod- 
ification reaction remain unknown. In this study, we 
mainly used small-angle X-ray scattering to charac- 
terize these conformational changes in solution and 
to unravel the mode of interaction between MnmE, 
MnmG and tRNA. In the nucleotide-free state MnmE 
and MnmG form an unanticipated asymmetric a.2p2 
complex. Unexpectedly, GTP binding promotes fur- 
ther oligomerization of the MnmEG complex leading 
to an a4p2 complex. The transition from the a2p2 to 
the a4p2 complex is fast, reversible and coupled to 
GTP binding and hydrolysis. We propose a model in 
which the nucleotide-induced changes in conforma- 
tion and oligomerization of MnmEG form an integral 
part of the tRNA modification reaction cycle. 

INTRODUCTION 

Transfer ribonucleic acid (tRNA) molecules contain a vast 
number of modified nucleotides. To date, over 90 of these 
modifications are known ranging from simple methyla- 
tions to complex hypermodifications (1,2). Those modifica- 



tions play structural or functional roles contributing to (i) 
the proper fold and stability of tRNA, (ii) proper codon- 
anticodon interaction at the decoding center of the ribo- 
some and (iii) tRNA recognition by the cognate aminoacyl- 
transferase (3). One of the main modification sites of tRNA 
is position 34, the so-called wobble position, that directly 
interacts with the third nucleotide of the messenger RNA 
(mRNA) codon. Considering their role in translation effi- 
ciency and fidelity, wobble modifications probably belong 
to the minimal set of tRNA modifications used in ancestral 
organisms (4). 

In bacteria, the proteins MnmE and MnmG 
form an enzyme complex (MnmEG) that is im- 
plicated in the modification of the wobble uri- 
dine in tRNALy^^nm5s2UUU, tRNA°'"mnm5s2UUC, 
tRNA "cmnm5s2UUG, tRNA '^"cmnmSUmAA, 

tRNA^'Sn^nmsucu and tRNAGiy^nmsucc (5,6). Except 
for the latter, all these tRNAs are reading A- and G-ending 
codons in split codon boxes (7,8). Depending on the sub- 
strate that is being used, the MnmEG complex first forms 
either 5-carboxymethylaminomethyluridine (cmnm^U- 
using glycine as substrate) or 5-aminomethyluridine 
(nm^U-using ammonium as substrate) (9). In a later step 
the bifunctional enzyme MnmC can convert these products 
to 5-methylaminomethyluridine (mnm^U), and finally the 
sulfur adding enzyme MnmA, in collaboration with a 
number of other proteins, will add a sulfur at position 2 of 
certain tRNAs, leading to mnm^s^U (10,11). tRNA^^^AA 
is an exception, as it does not get modified by either MnmC 
or MnmA, but it does get modified by TrmL, leading 
to the formation of 5-carboxymethlyaminomethyl-2'-0- 
methyluridine (cmnm^Um) (6). 

In eukaryotes, the orthologs of MnmE and MnmG 
are targeted to mitochondria and modify mitochondrial 
tRNAs (12). Interestingly, in human mitochondria, these 
orthologs (called GTPBP3 and MTOl, respectively) in- 
corporate during the modification reaction a taurine 
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molecule instead of glycine, leading to 5-taurinomethyl- 
uridine (Tm^U) (13,14). In bacteria, MnmE and especially 
MnmG have been identified as important regulators and 
determinants of bacterial virulence (15,16). In human, on 
the other hand, mutations of these enzymes are involved in 
severe mitochondrial myopathies (MELAS and MERRF) 
as well as in non-syndromic deafness (17), and the for- 
mer two diseases are known to be related to deficiencies in 
Tm^U tRNA modification (18). Moreover, it has been re- 
cently shown that mutations in MTOl cause hypertrophic 
cardiomyopathy and lactic acidosis (19). 

MnmE (formerly known as TrmE) is a homodimeric pro- 
tein of about 50-kDa subunits, where each subunit consists 
of an N-terminal domain, a helical domain and a G do- 
main that is inserted within the helical domain. The N- 
terminal domain is involved in homodimerization and is 
responsible for the binding of a tetrahydrofolate (THE) 
derivative. This THE derivative has been proposed to be a 
5 , 1 0-methylene-THE (MTHE) that serves as the one carbon 
donor for the C5 methylene moiety incorporated in uracil 
(9). MnmE belongs to the family of G proteins activated by 
nucleotide-dependent dimerization (GAD) (20,21). Com- 
pared to canonical small G proteins from the Ras family, 
GADs such as MnmE show a fast dimerization-dependent 
GTP hydrolysis rate combined with a low affinity for 
guanosine-5'-diphosphate (GDP) (22). This makes them in- 
dependent of guanine nucleotide exchange factors (GEEs) 
or GTPase activating proteins (GAPs) to cycle between a 
GTP-bound 'on state' and a GDP-bound 'off state' (20,21). 
Crystal structures of MnmE have shown that in the GDP- 
bound state the nucleotide-binding sites of the G domains 
are facing each other, but do not interact (22,23). How- 
ever, crystal structures of the isolated G domains show 
that these dimerize upon binding of stable GTP analogues 
or transition state analogues in a K"^ -dependent manner 
(24). The occurrence of such nucleotide-induced conforma- 
tional changes in the G domains is further supported in the 
context of the full-length protein by electron paramagnetic 
resonance (EPR) measurements (25). The dimerization of 
the G domains with concomitant reorganization of switch 
loops and catalytic machinery, together with the binding of 
a K"^ ion in the active site to stabilize the negative charges in 
the transition state, finally leads to hydrolysis of GTP with 
rate -limiting dissociation of the G domains (26). Previous 
experiments have shown that, in contrast to the classical G 
protein paradigm, an active GTP hydrolysis is required for 
the tRNA modification to occur in vitro and in vivo (23). 

MnmG (formerly known as GidA) is a homodimeric pro- 
tein consisting of 70-kDa subunits that bind flavin ade- 
nine dinucleotide (EAD) and nicotinamide adenine di- 
nucleotide. Crystal structures of MnmG from several or- 
ganisms have shown that each MnmG subunit consists of 
an EAD-binding domain, an insertion domain and a he- 
lical domain (27-29). Previous experiments indicated that 
MnmG is mainly responsible for tRNA binding within the 
MnmEG complex (28). Furthermore, two conserved cys- 
teine residues in the vicinity of the active site were identified 
to be crucial for the tRNA modification reaction and are 
proposed to play a catalytical role (28). 

The currently proposed mechanism of the MnmEG- 
catalyzed tRNA modification reaction implies that the 



activities of MnmE and MnmG are interdependent and 
that co-factors from both proteins (MTHE and EAD) are 
needed simultaneously in the tRNA modification reaction 
(6,9). We have previously proposed that the large confor- 
mational changes of the G domains of MnmE upon GTP 
binding and hydrolysis are relayed throughout MnmE and 
MnmG and orchestrate the tRNA modification reaction 
(23,25). However, awaiting structural information of the 
MnmEG complex in different nucleotide-bound states, the 
nature and relevance of these conformational changes re- 
main poorly understood. 

Here we used small-angle X-ray scattering (SAXS) to un- 
ravel the mode of interaction between MnmE and MnmG 
in the a2(32 complex (i.e. one MnmE dimer bound to one 
MnmG dimer). Surprisingly, MnmE and MnmG interact 
in an asymmetric manner which is distinct from the ear- 
lier proposed model. This model is validated using biophys- 
ical measurements. Size exclusion chromatography (SEC) 
experiments as well as SEC coupled to multiangle light 
scattering (SEC-MALS) and SAXS indicate that MnmE 
and MnmG form a higher oligomeric state when GTP is 
bound (a4(32, i.e. two MnmE dimers binding to one MnmG 
dimer), and this oligomerization appears to be reversible 
upon GTP hydrolysis. We propose a model in which the 
nucleotide-induced changes in conformation and oligomer- 
ization of MnmEG form an integral part of the tRNA mod- 
ification reaction cycle. 

MATERIALS AND METHODS 

Protein expression and purification 

The open reading frame coding for MnmE from Escherichia 
coli was cloned in a pET20 vector and the open reading 
frames coding for MnmG from E. coli and Aquifex aeoli- 
cus were cloned in a pET14b and a pET28a vector, respec- 
tively, as described previously (27). Point variants of MnmE 
and MnmG were prepared via the QuikChange mutagen- 
esis protocol from Stratagene. All proteins were expressed 
in Rosetta(DE3)pLysS and purified as described earlier 
(27). Eor preparation of the MnmE-MnmG a2(32 complex, 
equimolar amounts of MnmE and MnmG were rapidly 
mixed in the absence of nucleotides. Eor preparation of the 
MnmE-MnmG-MnmE a4(32 complex, MnmE was incu- 
bated with 1-mM GDP-AlEx (1-mM GDP, 1-mM AICI3 
and 10-mM NaE) and subsequently MnmG was added ei- 
ther in equimolar amount or in a 2 MnmE: 1 MnmG stoi- 
chiometry. 

Generation and purification of MnmG-specific nanobodies 

A llama was immunized over a period of 6 weeks with in to- 
tal 2 mg of purified MnmEG complex. The immunization, 
library construction and selection have been performed fol- 
lowing standard procedures (30) with minor modifications: 
total RNA was extracted from the peripheral blood lym- 
phocytes (31) and 50 ixg of total RNA was used to pre- 
pare complementary deoxyribonucleic acid (cDNA) using 
Superscript II (Invitrogen) and dN6 random primers, ac- 
cording to the manufacture's instruction. MnmEG-specific 
phage was eluted from the MnmEG-coated wells with 100- 
mM triethylamine (pH 10) and used to infect fresh TGI 
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cells. The MnmEG-coated wells were washed once with 
Tiis-HCl (pH 6.8), several times with phosphate buffered 
saline, and freshly grown TGI cells were added to the 
wells to recover the non-eluted phage. Sequence analysis 
revealed nine different nanobody families. Enzyme-linked 
immunosorbent assay confirmed that most of these, includ- 
ing Nb_MnmG_l, recognize MnmG Finally, all selected 
nanobody genes were cloned in a pHEN6 vector for expres- 
sion with a histidine-tag in E. coli and nanobodies were pu- 
rified via standard procedures (30). 

Preparation of tRNA 

The genes coding for E. coli tRNA^y^(UUU) and A. aeoli- 
cus tRNA^y**(UUU) were cloned in a pAsk-Iba vector and 
the tDNA including the T7 promotor site was amplified via 
polymerase chain reaction (PGR). tRNA was transcribed 
from the PGR product using a T7 RNA polymerase kit 
(Promega) via run off transcription according to the man- 
ufacturer's instructions. Unincorporated nucleotides were 
removed using MicroSpin G-25 columns (GE Healthcare) 
and tRNA was purified via precipitation using 0.1 volumes 
of 3-M sodium acetate and 2.5 volumes of 95% ethanol. To 
ensure folding of the tRNA to its global energy minimum, 
the tRNA transcript was heated to 368 K and slowly cooled 
to 288 K in a buffer consisting of 10-mM MgCli, 100-mM 
KCl and 20-mM Tris-pH 7.5. 

SAXS data collection and analysis 

SAXS experiments were performed at the beamlines P12 
and X33 of the EMBL in Hamburg (Germany), at the 
Swing beamline in Soleil (Paris, France) and at the ID 14- 

3 BioSAXS beamline at the ESRF in Grenoble (France). 
In the cases where a batch setup was used, scattering inten- 
sities were collected at protein concentrations of 8 mg/ml, 

4 mg/ml, 2 mg/ml and 1 mg/ml. In the cases where an 
online high-performance liquid chromatography (HPLC) 
setup was used, scattering intensities were collected after 
injection of 50 |xl of 8-mg/ml protein or protein complex 
on a KW404 gel filtration column (Shodex). This allowed 
us to collect online SAXS data on the elution peak of in- 
terest. In case an online stopped-flow setup was used, 32.5- 
|jlM MnmE and 32.5-jjlM MnmG were rapidly mixed in the 
measuring capillary in a 1 : 1 ratio, using a Bio-Logic SFM- 
300 stopped-flow apparatus coupled to the P12 beamline 
(EMBL, DESY, Hamburg) and data were collected at dif- 
ferent time intervals after mixing. As controls, either 32.5- 
|jlM MnmE or 32.5-|xM MnmG was mixed with buffer. All 
SAXS measurements were performed in 20-mM Hepes pH 
7.5, 150-mM KCl, 5-mM MgCh and 2-mM dithiothreitol 
(DTT), either with or without added nucleotides depend- 
ing on the experiment. Data were processed using PRIMUS 
(32). The radius of gyration {R^) was evaluated using the 
Guinier approximation (3 3 ) and also from the entire scatter- 
ing curve using Porod's law (34), and the distance distribu- 
tion function p{r) was calculated using gnom (35). Guinier 
plots of all SAXS data support monodispersity of the an- 
alyzed samples (Supplementary Figure SI). The measured 
scatter curves were compared with the theoretical scatter- 
ing curves of the protein models using Crysol (36). Ab ini- 



tio shape reconstructions were performed using 10-20 inde- 
pendent runs of Dammif (37) followed by Damaver (38) to 
create the final ah initio shape. OLIGOMER (39) was used 
to evaluated the percentage of open and closed MnmE in 
the GppNHp-bound state. ScAtter was used to calculated 
the molecular weights of the proteins and for the calculation 
of the fitting parameter i?sAS, the small-angle scattering in- 
variant Vc and the parameter (40). The summary of all 
statistics is given in Supplementary Table SI. 

Protein modeling, docking and refinement 

A model of E. coli MnmE was prepared by homol- 
ogy modeling starting from Chlorobium tepidum MnmE 
[PDB code 3GEE (25)] using MODELLER (41) after per- 
forming sequence alignments with Kalign (42). A model 
of Nb_MnmG_l was prepared starting from a canonical 
nanobody directed against carbonic anhydrase [PDB code 
1F2X (43)] using MODELLER. The complenentarity de- 
terimining regions (CDRs) were refined using GROMACS 
(44). Docking of MnmE to MnmG, tRNA to MnmG and 
Nb_MnmG_l to MnmG was performed using PatchDock 
(45,46). Hereby, the two partners are docked based on shape 
matching algorithms and the docking solutions are ranked 
according to a geometric shape complementarity score. 
Subsequently, the theoretical scattering curves of all 413 
docking models were compared to the experimental SAXS 
curve to filter out the best fitting model (see Supplementary 
material and methods for a detailed description of the selec- 
tion of the MnmEG a2p2 complex). Missing residues not 
seen in the crystal structure were added using MODELLER 
(41) and the N-terminal tag was added in Coot (47). Models 
were refined using molecular dynamics simulations in water 
using the OPLS force field within GROMACS (44). All sim- 
ulations were done at a constant temperature of 300 K and 
a constant pressure of 1 atm over a simulation time of 1 ns. 
Final figures were created in PyMOL (www.pymol.org). 

Isothermal titration calorimetry 

Isothermal titration calorimetry (ITC) experiments were 
performed at 278 K using the MicroCal iTC200 sys- 
tem (GE Healthcare) in a buffer consisting of 20-mM 
Hepes pH 7.5, 150-mM NaCl, 5-mM MgCli and 2-mM 
(3-mercaptoethanol (BME). For binding measurements of 
wild-type MnmE and MnmG, a protein concentration in 
the cell and syringe of, respectively, 15 |jlM and 300 |xM 
was used, in the presence of 1-mM FAD. Titrations were 
performed in both directions with either MnmE or MnmG 
in the measuring cell. The binding constants of variants of 
MnmE and MnmG were measured by titration of 700-|xM 
MnmE (in the syringe) into 50-|xM MnmG (in the cell) and 
vice versa, in the presence of 1-mM FAD. The resulting data 
were fitted to a one site binding model using the Origin soft- 
ware accompanying the ITC instrument. 

Analytical high-resolution size exclusion chromatography 

The formation of complexes between MnmE and MnmG, 
as well as with tRNA, was followed using high-resolution 
analytical gel filtration experiments. Proteins and tRNA 
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were mixed and, after 20-min incubation, 25 |xl was in- 
jected on a KW803 gel filtration column (Shodex) cou- 
pled to an HPLC system (Waters). The separation was per- 
formed at a flow rate of 0.8 ml/min in a buffer consisting 
of 50-mM Hepes pH 7.5, 150-mM KCl, 5-mM MgCh and 
5-mM BME. For the complex formation in the presence of 
GDP-AlFx, MnmE was incubated with 1-mM GDP-AlFx 
(1-mM GDP, 1-mM AICI3 and 10-mM NaF) for 30 min and 
then mixed with MnmG. 

GTP hydrolysis assay 

Hydrolysis of GTP by the MnmE(E282A)-MnmG com- 
plex was followed over time by separation of the nu- 
cleotides (GTP and GDP) using a C18 reversed phase col- 
umn attached to an HPLC system (Waters). 66.7-|xM E. 
coli MnmE(E282A)-MnmG complex was mixed with 100- 
|jlM GTP in a buffer consisting of 20-mM Hepes pH 7.5, 
150-mM KCl, 5-mM MgClj and 2-mM DTT At several 
time points, between 0 min and 24 h, 20-|xl samples were 
taken and the reaction was stopped by flash freezing in 
liquid nitrogen. After boiling and centrifugation, the su- 
pernatant was loaded on a CI 8 reversed phase column 
(Jupiter, 25 cm x 4.6 mm) and the nucleotides were eluted 
using a buffer containing 100-mM KH2PO4 pH 6.4, 10-mM 
tetrabutyl ammonium bromide and 7.5% acetonitrile. Nu- 
cleotides were detected using absorbance at 254 nm. The 
GTP concentration in the sample was calculated from the 
peak area using a standard curve derived from known GTP 
concentrations. At each time point, also a second aliquot 
was taken, which was immediately analyzed using analyt- 
ical high-resolution size exclusion chromatography (as de- 
scribed above), to assay complex formation between MnmE 
and MnmG. 

RESULTS 

Dimerization of the G domains of MnmE coincides with 
movements of the helical domains 

X-ray crystal structures of MnmE have until now been re- 
ported from C. tepidum, Nostoc and Thermotoga maritima 
(25). These structures show MnmE in the so called 'open' 
state, with the G domains of the homodimer facing — but 
not contacting — each other. However, crystal structures of 
the isolated G domains oiE. coli MnmE bound to the tran- 
sition state homologue GDP-AlFx showed that the G do- 
mains dimerize in the transition state (24). Supported by 
EPR studies and biochemical studies (25), this led to a 
model for MnmE action in which the protein cycles between 
a GDP-bound 'open' state and a GTP-bound 'closed' state. 
Unfortunately, until this date, no structures have been re- 
ported of a full-length MnmE in this closed state. In a first 
stage therefore, we set out to investigate the conformational 
behavior of E. coli MnmE in solution using SAXS measure- 
ments. Since currently no crystal structures of full-length E. 
coli MnmE are available, a homology model of this protein 
in the open state was generated starting from the existing 
crystal structure of C. tepidum MnmE (PDB code 3GEE) 
(25), with which it shares 36% sequence identity (Figure 
lA). 



SAXS data of the 'open' E. coli MnmE were collected 
in batch mode at the X33 beamline (Hamburg, Germany) 
in the absence of nucleotides. Analysis with Crysol and 
ScAtter shows that the scattering curve of E. coli MnmE 
agrees very well with the homology model (chi-: 0.82; -Rsas^ 
0.0002) validating the crystal structures of the open form 
of MnmE in solution (Figure IC). Also the theoretical Rg 
(radius of gyration) of the model agrees well with the experi- 
mental Rg (-Rgmodei: 38.1 A; i?gexp: 38.6 A ± 0.9 A). Further- 
more, the calculation of the molecular weight using either 
the SAXS data (MWsaxs) or SEC-MALS data (MWmals) 
(Supplementary Figure S2 and Supplementary Table S2) 
supports that E. coli MnmE is a dimer in solution (MWcaic^ 
103 kDa; MWsaxs: 104 kDa; MWmals: 100 kDa). Subse- 
quently, an ah initio molecular envelope was created from 
the experimental SAXS data using Dammif This envelope 
overlays perfectly with the homology model supporting that 
the G domains of MnmE are separated in the open state 
(Figure lA). 

To gain insights into the structure of the closed state of 
E. coli MnmE, SAXS data were collected in the presence 
of 1-mM GDP-AlFx at the X33 beamline (Hamburg, Ger- 
many) in batch mode. Clearly distinct scattering curves are 
obtained for MnmE in the nucleotide -free (open) state and 
the GDP-AlFx-bound (closed) state (Figure IC), indicative 
of considerable conformational changes in solution. Since 
no crystal structures were available for the closed state, we 
constructed a model, whereby the G domains were arranged 
as found in the crystal structure of the isolated G domains 
in the closed state (24), but leaving the other domains of 
MnmE unchanged (see Figure IB). This model fits the ex- 
perimental data reasonably well (chi-: 2.4; -Rsas: 0.0001) 
and is clearly a better approximation of the structure in so- 
lution than the open form of MnmE, thus further support- 
ing occurrence of the open-close transition of MnmE (Sup- 
plementary Figure S3). The Rg of the model agrees well with 

the experimental Rg {Rg^oieV- 37.0A; i?gexp: 37.2 A ± 0.9 A) 
and the calculation of the molecular weight supports that 
E. coli MnmE is also dimeric in the closed form (MWcaic: 
103 kDa; MWsaxs: 121 kDa; MWmals: 97 kDa). Based on 
the experimental SAXS data an ab initio shape was calcu- 
lated using Dammif (Figure IB). This molecular shape con- 
firms the dimerization of the G domains but suggests that 
an additional 'upward' movement of the helical domains of 
MnmE is taking place concomitant with the closing of the 
G domains (indicated by arrows in Figure IB). 

EPR studies have previously shown that, while MnmE 
mainly exists in the open state if GDP is bound and in the 
closed state if GDP-AlFx is bound, an equilibrium between 
both states exists in the presence of GppNHp, a ground- 
state analog of GTP (25). Correspondingly, an experimental 
SAXS curve of MnmE in the presence of 1 mM of GppNHp 
could neither be appropriately fitted to the model of the 
open state nor to the closed state. Therefore OLIGOMER 
was used to estimate the fractions of the open and closed 
forms of MnmE in the GppNHp-bound state (Figure ID). 
This approach led to an optimal description of the SAXS 
data (chi^: 1.8; 7?sas: 0.0002) using a ratio of 88% MnmE in 
the closed state and 12% in the open state. This fraction of 
MnmE in the closed form is higher than that found via EPR 
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Figure 1. Model of £. coli MiimE in the (A) open and (B) closed states superposed on their ah initio shape as obtained by Dammif. The closure of the G 
domains probably also induces an 'upward' shift in the helical domains as indicated by the arrows. (C) Superposition of experimental and theoretical scatter 
curves using Crysol. The experimental scattering curves of the open (green) and the closed (blue) state of MnmE agree well with the theoretical scattering 
curves (red) obtained from the models. The residuals of the fitting are shown below the experimental curves in the corresponding color. (D) Fitting of the 
experimental scatter curve of MnmE bound to GppNHp to a mixture of MnmE in the open and closed states using OLIGOMER. A fraction of 88% of 
MnmE in the closed state and 12% of MnmE in the open state is obtained. The residuals of the fits are depicted below. 



(about 30%), which might be due to the different conditions 
and temperatures that were used in both experiments. 

The solution structure of MnmG agrees well with the crystal 
structure 

Currently crystal structures of MnmG have been reported 
from E. coli, A. aeolicus and C. tepidum (27-29). To confirm 
whether these crystal structures correspond to the protein 
structure in solution, SAXS data were collected. Unfortu- 
nately, slight aggregation of E. coli MnmG in the capillary 
hampered accurate interpretation of the scattering data. 
Therefore we turned to the homologous MnmG from A. ae- 
olicus (51% sequence identity to E. coli MnmG) (PDB code 
2ZXI) (28). Surprisingly, SEC-MALS data suggest that 



MnmG from A. aeolicus is mainly present as a monomer 
in solution, with only a minor dimer peak at the concen- 
trations used [Supplementary Figure S2(F)]. However, since 
E. coli MnmG is clearly dimeric at all concentrations used 
[Supplementary Figure S2(D)], we focused our SAXS anal- 
ysis on the dimer peak of A. aeolicus MnmG. Therefore, 
SAXS data were collected at the SWING beamhne (Soleil, 
Paris) using an HPLC setup, in which purified MnmG is ap- 
plied on a gel filtration column and the SAXS spectra were 
recorded after elution of the protein from the gel filtration 
column. The theoretical scattering curve of the A. aeolicus 
MnmG dimer (as observed in the crystal structures) shows 
a good fit to the experimental SAXS data (chi^: 1.9; i?sAs: 
0.03) (Supplementary Figure S4). The obtained Rg values 
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(-Rgmodei: 40.2 A; 7?gexp: 42.3 A ± 1.5 A) and the calculated 
molecular weights further support that A. aeolicus can also 
form a dimer in solution (MWcaic: 150 kDa; MWsaxs: 168 
kDa; MWmals: 130 kDa). An ab initio molecular model 
was subsequently created from the data using Dammif This 
envelope can be nicely superimposed on the crystal struc- 
ture (Supplementary Figure S4). 

The MnmG dimer seems to bind a single tRNA molecule 

Within the context of the tRNA modification reaction, 
it has been shown that MnmG is mainly responsible for 
tRNA binding and can bind tRNA in the absence of MnmE 
(28). To prepare an MnmG-tRNA complex, we mixed 
A. aeolicus MnmG with a 1.2-molar excess of A. aeolicus 
tRNA^^'s^UUU) and purified the complex via gel filtration. 
SAXS data of the fractions containing the MnmG-tRNA 
complex were afterward measured at the SWING beam- 
line (Soleil, Paris) using an HPLC setup. In order to create 
a model of the MnmG-tRNA complex, first docking was 
performed using PatchDock. The MnmG-tRNA docking 
models were subsequently selected against the experimen- 
tal SAXS data based on the chi^ value obtained by Crysol. 
This procedure led to two equally well scoring models that 
cannot be distinguished by SAXS (see Figure 2 and Sup- 
plementary Figure S5 for a comparison of the models). For 
both models, however, the data clearly suggest that only 
one subunit of MnmG is bound to a tRNA molecule, while 
the second subunit remains vacant. This is established by 
the good fit of the experimental data to a model with only 
one bound tRNA (chi-: 2.8; -Rsas: 0.006), in comparison 
to a similar model with two bound tRNAs (chi-: 22; i?sAs: 
0.03) (Figure 2B and Supplementary Figure S6). Addition- 
ally, the ab initio shape created based on the experimental 
data with Dammif shows extra features on only one side 
of the MnmG molecule, in agreement with the presence of 
only one tRNA (Figure 2A). However, it might be possi- 
ble that we did not saturate the binding sites of MnmG and 
that both tRNA-binding sites might still be occupied when 
using higher amounts of tRNA. In addition, we therefore 
also performed SEC-MALS experiments on the A. aeolicus 
tRNA and the A. aeolicus MnmG-tRNA complex. To this 
end, we mixed a 5-fold molar excess of tRNA with MnmG 
before injection on the size exclusion column. The molecu- 
lar weight obtained for the A. aeolicus MnmG-tRNA com- 
plex is 151 kDa [see Supplementary Figure S2(H)], which 
is somewhat lower than the calculated molecular weight 
of one tRNA molecule binding to an MnmG dimer (174 
kDa), but which corresponds to the sum of the experimental 
molecular weight from SEC-MALS of MnmG (130 kDa) 
and tRNA (25 kDa) [see Supplementary Figure S2(F) and 
(G)]. 

Despite the fact that we are unable to discriminate be- 
tween the two best scoring MnmG-tRNA models (Supple- 
mentary Figure S5) solely based on the scattering data, we 
clearly prefer the model shown in Figure 2C based on mech- 
anistic arguments. In this model, the anticodon stem-loop 
(ASL) of the tRNA, containing the substrate wobble uri- 
dine, is located in the FAD-binding pocket of MnmG close 
to the co-factor FAD and to the catalytic residues C48 and 
C248 {A. aeolicus numbering), of which the latter is sit- 



uated at the end of a flexible loop [compare Supplemen- 
tary Figure S5(C) and (D)] (28). This MnmG-tRNA model 
is further supported by a previous site-directed mutagene- 
sis study that was used to map the tRNA-binding region 
of A. aeolicus MnmG (28). The tRNA-interacting residues 
that were identified in this way (e.g. Asn49, Arg97, Arg282, 
Lys290, Lys293, Arg439, Arg443, Tyr545 and Arg548; A. 
aeolicus numbering) coincide with the MnmG-tRNA inter- 
action surface in our model (see Figure 2C). 

MnmE and MnmG form an asymmetric, L-shaped complex 

As already shown previously, MnmE and MnmG can 
form a tetrameric complex with a2(32 stoichiometry (i.e. 1- 
MnmE and 1-MnmG homodimer) in the GDP-bound state 
(23). Complex formation can be monitored by size exclusion 
chromatography, where a 1 : 1 mixture of MnmE and MnmG 
at 100 |xM elutes at a smaller elution volume than the in- 
dividual proteins, corresponding to a molecular weight of 
260 kDa, close to the expected molecular weight (243 kDa) 
for the MnmE-MnmG a2(32 complex (Figure 3A). How- 
ever, upon mixing of MnmE and MnmG at lower concen- 
trations a gradual shift in the elution volume is observed 
from a volume corresponding to the molecular weight of the 
a2p2 complex to the elution volume corresponding to the 
individual proteins. This is typical for a low affinity complex 
with a very fast binding and dissociation rate. Moreover, at 
very high concentrations of MnmE and MnmG an even fur- 
ther decrease in elution volume, corresponding to a higher 
molecular weight complex, can be observed. 

To get insights into the overall architecture of the a2p2 
complex E. coli MnmE was mixed equimolar with E. coli 
MnmG and the resulting complex was purified via gel fil- 
tration. SAXS data were collected on this sample at the 
SWING beamline (Soleil, France) using an HPLC setup. 

A model of the MnmEG complex was generated by dock- 
ing of MnmG to MnmE. Hereto, we started from MnmE 
with its G domains in the closed form, since both EPR 
studies and stopped-flow fluorescence studies indicated that 
MnmG binding at least partially induces G domain dimer- 
ization in MnmE (23, 48). The resulting docking models 
were validated and ranked using experimental SAXS data, 
based on the chi^ value obtained by Crysol and on agree- 
ment with the ab initio shape as obtained by Dammif, and 
finally on visual inspection of the models in terms of func- 
tionality (proximity of active sites, vacancy of the tRNA- 
binding site on MnmG — see Supplementary material and 
methods for a detailed description). This led to the model 
shown in Figure 3C and D. Considering the size of the com- 
plex and the fact that we cannot exclude the occurrence 
of conformational changes in both proteins upon complex 
formation, a good fit of the theoretical and experimental 
scatter curves is obtained (chi^: 3.1; i?sAs = 0.0005) (Fig- 
ure 3B). The 1:1 stoichiometry of the model is supported by 
the radius of gyration (i^gmodei^ 51.2 A; i?gexp: 52.3 A ± 1.4 

A) and the molecular weight (MWcak: 243 kDa; MWsaxs: 
207 kDa) and is also further confirmed by the molecular 
weight obtained by SEC-MALS measurements (MWmals: 
213 kDa; Supplementary Figure S2). The fact that the ex- 
perimental molecular weights are somewhat smaller than 
the calculated molecular weight probably indicates that at 
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Figure 2. (A) Docking model of the A. aeolicus MnmG-tRNA complex superposed on the ah initio shape as obtained by Dammif. (B) The fit of the 
theoretical (red) to the experimental (blue) scattering curve supports the model. The residuals of the fits are also shown. (C) Details of the preferred model 
of the MnmG-tRNA interaction (compare to an alternative model in Supplementary Figure S5) and (D) a close-up view of the (FAD-binding) active site. 
This model is supported by the proximity of known tRNA-interacting residues (red) (28) and the catalytic cysteines (C48 and C284; yellow) and by the 
location of the FAD co-factor (shown as a stick model) with respect to the tRNA wobble nucleotide (green). 



the concentration used, the complex formation is not com- 
pletely saturated. 

In this model, one MnmE dimer is binding via the N- 
terminal domain and the helical domain of one subunit to 
the C-terminal domain (mainly helix otll and a 19) of one 
subunit of the MnmG dimer in a non-symmetric manner 
(Figure 3C). Interestingly, this model leaves one protein- 
binding site vacant on MnmE as well as one on MnmG, 
allowing for further oligomerization. The fact that such 
oligomerization is not observed below 200 |xM (in the ab- 
sence of nucleotides, see further) indicates a form of neg- 
ative cooperativity between the protomers of both MnmE 
and MnmG. A similar negative cooperativity was also ob- 
served between the subunits of MnmG with regard to tRNA 
binding. However, contrary to the previously proposed 
models for the MnmE-MnmG interaction (27), in our cur- 
rent model the 5,10-methylenetetrahydrofolate-binding site 
of MnmE and the FAD-binding site of MnmG are oriented 
toward each other. Such an arrangement would be required 



to be conform with the proposed mechanism of tRNA mod- 
ification by MnmEG (6, 9), which assumes an intensive col- 
laboration between both active sites. Finally it should be 
noted that the MnmEG interaction model that was previ- 
ously proposed (27) does clearly not fit with the experimen- 
tal SAXS data (chi^: 30.3; i?sAs = 0.01) (Supplementary 
Figure S7). 

In order to validate our model we used ITC measure- 
ments. Wild-type MnmE binds MnmG with an apparent 
dissociation constant {Kx)) of 3 |jlM and 7 jjlM, depending 
on the direction of the titration (Figure 4). To our surprise, 
we observed that the stoichiometry of the complex con- 
sistently depends on the direction of the titration. Adding 
MnmG (in the syringe) to MnmE (in the cell) leads to a stoi- 
chiometry of 1 , which is consistent with a 1 : 1 complex (23). 
However adding MnmE to MnmG leads to a stoichiome- 
try of 0.5. This observation was repeatable with different 
batches of protein, but we do not have a fully adequate ex- 
planation so far. However, such a phenomenon has been ob- 
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Figure 3. (A) Gel filtration experiment (KW-803, Shodex) showing the complex formation of the MnmE dimer (blue) and the MnmG dimer (red), forming 
the c(2p2 MnmEG complex (black). Twenty five microliters of 100-jjlM solutions were injected on the column. (B) The theoretical scattering curve (red) 
of the MnmEG docking model agrees well with the experimental scattering curve (blue). The residuals of the fits are also shown. (C) and (D): Docking 
model of the a2p2 MnmEG complex superposed on the ab initio shape as obtained by Dammif (gray). 



served before in the case of systems capable of forming long 
alternating chains of homodimers, like ccdA/ccdB (49). In 
that specific case, it was indicative of kinetic limitations that 
prevent reaching equilibrium in at least one direction of the 
titration. Here, a full interpretation only became possible 
after high resolution structures of the ccdA/ccdB protein 
complexes were solved. 

Subsequently we measured the effect on the binding affin- 
ity of point mutations of conserved residues in the proposed 
binding interface (Figure 4 and Supplementary Figure S8). 
It should be stressed here that the SAXS-based modeling 
only provides low-resolution shapes of the complexes be- 
tween MnmE and MnmG and thus any detailed informa- 
tion concerning individual atomic interactions is missing. 
As pointed out before we can also not exclude the occur- 
rence of small conformational changes in the individual 
proteins upon complex formation. To more easily pinpoint 
the interaction surfaces by mutagenesis we therefore chose 
to construct multiple mutants, hence covering a whole patch 
on the proposed interaction surface. To exclude the possi- 
bility that these mutations have an influence on the over- 
all structure and stability of the proteins, rather than dis- 
rupting the actual binding interface, we performed thermal 
unfolding experiments (Supplementary Figure S9). All mu- 
tants show a melting temperature that is comparable to the 
corresponding wild-type proteins, in the range of 47° C. As 
the ITC experiments were performed at 5°C, these data sup- 
port that any effect of the mutations on binding is not due to 
a destabilizing effect on the proteins. The strongest decrease 



in (apparent) affinity is observed for the mutants MnmE 
(N118A, K120A and D122A) and MnmG (Y551 A, R554A, 
E558A and R564A) that completely abolish binding at the 
concentrations used in our ITC experiment. Both patches 
of residues are located in the middle of the proposed Mn- 
mEG interface and are in very close proximity to each other. 
Mutations that are located at the periphery of the proposed 
interaction surface clearly decrease affinity to a lesser ex- 
tend. The mutant MnmE (Y169A, D175A and F176A), lo- 
cated near the outer end of the helical domain, shows a 2- 
4-fold decrease in affinity (depending of the direction of the 
titration). An additional mutation of Serl53 to alanine, lo- 
cated at the edge of the N-terminal and helical domain, de- 
creases the (apparent) affinity even further (up to 10-fold 
compared to wt MnmEG). Also the mutant MnmE (T434A 
and D440A), located between the N-terminal domain and 
the helical domain, decreased the (apparent) affinity toward 
MnmG (2-3-fold). Finally, in our model the C-terminal end 
of MnmG (amino acids 550-629) is also involved in the in- 
teraction with MnmE. It has previously been shown that 
deletion of this part of MnmG severely affects binding of 
MnmG to MnmE (29). Altogether these mutagenesis data 
support the MnmEG model proposed here. It is notewor- 
thy that all mutants of MnmE that showed an effect on the 
binding affinity also reversed the stoichiometry of titration 
compared to wild-type MnmEG. The underlying mecha- 
nism for this observation remains thus far obscure. 
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Figure 4. Validation of the MnmEG complex interface using ITC. (A) Residues that showed a decrease in affinity upon mutation to alanine are indicated 
as followed - red: MnmE (N118A, K120A and D122A); green: MnmE (Y169A, D175A and F176A) with additional S153A; yellow: MnmE (T434A and 
D440A); cyan: MnmG (Y551A, R554A, E558A and R564A). MnmE is shown with its surface in light and dark gray and MnmG is shown as a cartoon 
model in purple and blue. ITC data of the titration of (B) MnmE wt to MnmG wt and (C) MnmG wt to MnmE wt, showing the dependence of the binding 
stoichiometry on the direction of titration. (D) Summary of binding stoichiometry (n) and apparent values of MnmE/MnmG variants as obtained by 
ITC. n.b. = no binding detected. 



GDP-AlFx binding induces formation of an a4p2 complex 

To further assess the influence of the nucleotide state (GTP 
versus GDP) of MnmE on its binding to MnmG, the 
MnmEG complex was analyzed in the presence of GDP- 
AlFx. High-resolution size exclusion chromatography sug- 
gests the formation of an even larger MnmEG complex 
when MnmE is bound to GDP-AlFx, even at very low 
protein concentrations (5 |jlM) (Figure 5A). The apparent 
molecular weight of this complex according to gel filtration 
is 520 kDa, compared to 260 kDa determined for the a2(32 
complex. Interestingly, a single peak for this larger com- 



plex on gel filtration could only be observed if MnmE was 
used in 2-fold excess over MnmG. A 1 : 1 mixture of MnmE- 
GDP-AlFx and MnmG led to the formation of a peak for 
this high molecular weight complex in addition to several 
other peaks at lower molecular weight. 

To gain insights into the overall architecture and shape 
of this larger complex, SAXS data were collected on the 
MnmE-GDP-AlFx-MnmG complex. Therefore MnmE- 
GDP-AlFx was mixed in 2-fold molar excess with MnmG 
and SAXS data were collected at the SWING beamline 
(Soleil, France) in an HPLC setup. The molecular weight as 



Nucleic Acids Research, 2014, Vol. 42, No. 9 5987 




Figure 5. (A) MnmE and MnmG form a a4p2 complex in the presence of GDP-AlFx. A single species on gel filtration is only obtained if MnmE is used 
in 2:1 excess over MnmG (orange) while a mixture of species is obtained in a 1:1 stoichiometry mixing (cyan). Twenty-five microliters of 100-(jlM solutions 
were injected on the column. (B) The theoretical scattering curve of the a4p2 model (red) agrees with the experimental scattering curve (blue). The residuals 
of the fits are also shown. (C) and (D): Model of the 0(4(32 MnmEG complex superposed on the ah initio shape as obtained by Dammif (gray). 



determined by SAXS was 338 kDa, which is in agreement 
with the molecular weight estimate from SEC-MALS (326 
kDa) and corresponds to a 2:1 (two MnmE dimers and one 
MnmG dimer) complex (MWcaic^ 346 kDa). This molecu- 
lar weight estimate is smaller than the estimate obtained 
from gel filtration. However, the molecular weight estima- 
tions of SAXS and MALS are independent of the shape 
of the molecule, whereas in gel filtration the extended na- 
ture of the a4(32 complex (see following) probably leads to 
the aberrant elution. We created a tentative model with 2:1 
(MnmE:MnmG) stoichiometry by placing a second MnmE 
dimer on the vacant binding site of MnmG in our a2(32 
model (Figure 5C and D). The theoretical scattering curve 
of the resulting model is in very good agreement with the 
experimental scattering curve (chi^: 4.7; i?sAs: 0.0007), es- 
pecially considering that possible smaller conformational 
changes upon complex formation in MnmE and MnmG 
are not taken into account (Figure 5B). The model is fur- 
ther supported by the obtained from the model com- 
pared to the experimental Rg (-Rgmodei: 67.2 A; R^txp- 67.4 
A ± 1.3 A). Finally, the ah initio molecular envelope ob- 
tained by Dammif agrees very well with the a4p2 model 
and shows the elongated nature of this complex (Figure 5C 
and D). As a control and to rule out that this larger complex 
is formed by the addition of two MnmG dimers to the va- 
cant binding sites of one MnmE, we created a model with 
a 1:2 stoichiometry (MnmE:MnmG). The fitting statistics 



clearly indicate that this model is incorrect (chi^: 303; i?sAs: 
0.12; 7?gmodei: 57.9 A; i^gexp: 67.4 A ± 1.3 A) and the model 
is also not in agreement with the ab initio envelope obtained 
by Dammif (not shown). 

Attempts to obtain a homogeneous sample of the 
MnmEG-tRNA complex did not succeed. Therefore, to get 
insights into this MnmEG-tRNA complex we generated a 
tentative model based on the previous models (Supplemen- 
tary Figure SIO). Since two tRNA-binding interfaces in the 
a4p2 complex are available, we suspect that this complex 
will bind two tRNA molecules. Therefore we superposed 
the MnmG-tRNA model (Figure 2 A) onto the model of 
the MnmE-GDP-AlFx-MnmG complex (Figure 5C). This 
tentative model would create some overlap between tRNA 
and MnmE. However, these sterical clashes could be re- 
lieved by the conformational changes in the a-helical do- 
mains of MnmE that are anticipated upon GTP binding 
(see Figure IB). This low-resolution architecture places the 
wobble position of the tRNA in between the 5,10-MTHF- 
and FAD-binding pockets of MnmE and MnmG, as would 
be required for the modification reaction to occur. The pro- 
posed movement of the ct-helical domain of MnmE upon 
GTP binding might bring the 5,10-MTHF in even closer 
contact to the wobble uridine of tRNA. Interestingly it 
should be noted that the MTHF closest to the wobble uri- 
dine of tRN A is not coming from the MnmE protomer in- 
teracting via its a-helical domain with MnmG, but is com- 
ing from the second MnmE protomer. 
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Nb_MnmG_l binds MnmG and disrupts the MnmE-MnmG 
complex 

A camelid single domain antibody fragment (Nanobody), 
which we named NbJVInmG_l, was identified to bind 
E. coli MnmG. Interestingly, gel filtration experiments 
show that Nb_MnmG_l can disrupt the MnmEG complex, 
both in the a2p2 and a4|32 forms, when added in excess 
[Supplementary Figure Sll(A)]. This finding suggests that 
Nb_MnmG_l and MnmE are competing for the same bind- 
ing surface of MnmG and can hence be used to further val- 
idate our MnmEG model. 

An ensemble of MnmG-Nb_MnmG_l models was cre- 
ated by docking of an Nb_MnmG_l homology model to 
MnmG and the resulting models were selected based on 
the experimental SAXS data of the MnmG-Nb_MnmG_l 
complex. The resulting best model, with NbJVInmG_l 
bound to MnmG in a 1:1 stoichiometry, is in good agree- 
ment with the experimental scattering curve (chi-: 3.6; -Rsas: 
0.0004) and with the ah initio model created with Dammif 
[Supplementary Figure S 11(B) and (C)]. This is in agree- 
ment with the obtained R^s (-Rgmodei: 42.6 A; i?gexp: 43.3 A ± 
0.4 A) and the calculation of the molecular weight from the 
SAXS data, as well as from the SEC-MALS data (MWcaic: 
173 kDa; Mwsaxs: 184 kDa; MWmals: 170 kDa). In this 
model, Nb_MnmG_l is binding to the C-terminal helical 
domain of MnmG, on a position that indeed partially over- 
laps with the binding surface of MnmE in our MnmEG 
model. This finding hence further confirms the validity of 
our proposed MnmEG interaction model. 

Conversion of the a2|i2 and (x4p2 complexes is linked to 
GTP binding and hydrolysis 

In the previous experiments, we have observed that MnmE 
and MnmG mainly form an a4(32 complex upon binding of 
GDP-AlFx, while an a2(32 complex is formed upon bind- 
ing of GDP at the the concentrations used in the gel filtra- 
tion experiments. On the other hand, under physiological 
conditions MnmE will actively hydrolyze the bound GTP 
and it has been previously described that this GTP hydrol- 
ysis is required for tRNA modification (25). We were there- 
fore interested to see whether the conversion of the a4(32 
complex into the a2(32 complex is linked to GTP hydroly- 
sis. 

Since GTP hydrolysis by wild-type E. coli MnmE occurs 
on time scales that are too fast to follow via gel filtration ex- 
periments, a slow hydrolyzing mutant MnmE (E282A) was 
used (24). The MnmE(E282A)-MnmG complex was mixed 
with a small excess of GTP and both the GTP hydrolysis 
and the oligomeric state of the protein complex were fol- 
lowed and compared over time (Figure 6 A and B). At time 
point zero, the MnmE(E282A)-MnmG complex is mainly 
GTP bound and correspondingly, mainly the a4(32 complex 
is observed (although some a2(32 complex is observed due 
to the equimolar mixing of MnmE with MnmG) (Figure 
6B). After 24 h all the complex is converted to the a2(32 
form and at intermediate time points a gradual shift from 
the large to the smaller complex is observed. This conver- 
sion in oligomeric state nearly perfectly corresponds to the 
hydrolysis of GTP by MnmE (E282A), which shows that all 



GTP is hydrolyzed after 24 h with a 50% GTP conversion 
after roughly 2 h (Figure 6A). 

This clearly indicates that GTP hydrolysis and complex 
disassembly are directly linked. However, the latter also im- 
plies that for the wild-type MnmEG complex, formation of 
the a4(32 complex from the a2(32 complex will need to com- 
pete with the disassembly due to GTP hydrolysis [/ccat = 0.16 
s~' (23)]. In order to be physiologically relevant, the a4(32 
complex must hence form in a time scale faster than or com- 
parable to GTP hydrolysis. To gain insights into how fast the 
complexes assemble, a stopped-flow setup directly coupled 
to SAXS beamline P12 (EMBL, Hamburg) was used. In this 
setup, MnmE and MnmG were rapidly mixed in the pres- 
ence of GDP-AlFx and SAXS spectra were recorded at dif- 
ferent time points after mixing (50-ms exposures). As con- 
trol, mixing of either MnmE or MnmG with buffer resulted 
in time invariable Rg^ of 3.82 nm and 4.63 nm, respectively, 
very close to the values obtained in a static HPLC setup 
(3.87 nm and 4.36 nm, respectively). Mixing of MnmE with 
MnmG in the presence of GDP-AlFx led to the formation 
of a large complex with an Rg fluctuating around 6.5 nm, in 
agreement with the formation of the a4(32 complex (Figure 
6C). It must, however, be noted here that under the high pro- 
tein concentration that is required in this experiment (due to 
the very low exposure times), even in the absence of GDP- 
AlFx the equilibrium of the MnmEG complex seems to be 
shifted to the a4(32 complex (data not shown). In either 
case, formation of the a4(32 complex upon mixing of MnmE 
and MnmG occurs very fast and within the dead time (<20 
ms) of our measurements. This thus strongly suggests that 
the rate of formation of this complex can compete with GTP 
hydrolysis and concomitant disassembly of the complex and 
that the a4p2 complex would hence be populated and rele- 
vant under active GTP turnover conditions. 

DISCUSSION 

The survival of any organism depends on the efficient and 
correct translation of mRNA into protein. In this pro- 
cess, the proper recognition of the codon (especially in split 
codon boxes) is immediately linked to the modification at 
the wobble position of tRNA. Hence, defects in wobble 
modification are often correlated with diseases in humans 
(18). 

MnmE and MnmG are responsible for the cmnm^U34 
wobble modification and act in a concerted manner (23). 
Although crystal structures of the single proteins were be- 
ing unraveled starting from 2005, it has proven to be very 
difficult to understand how they collaborate in the tRNA 
modification reaction. It was known that the proteins could 
form an a2(32 complex (23), but no structures of the Mn- 
mEG complex nor of complexes with the tRNA substrates 
have been reported to date. 

Here, we used solution SAXS studies to reveal the in- 
teraction and interplay between MnmE and MnmG at 
low resolution. Based on our studies, we propose a model 
(Figure 7) for the different steps that occur during the 
tRNA modification reaction by MnmEG. Based on bio- 
chemical and EPR studies (24,25), it was previously sug- 
gested that MnmE undergoes a transition from an 'open' 
(GDP-bound) to a 'closed' (GTP-bound) conformation. 
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Figure 6. (A) GTP hydrolysis by the MnmE(E282A)-MnmG complex shows a similar time dependency to (B) the concomitant conversion of the 
MnmE(E282A)-MnmG complex from the a4p2 to the a2p2 form during GTP hydrolysis. (C) Analysis of MnmE-MnmG complex formation using 
stopped-flow SAXS. The evolution of the R^^ in function of time upon mixing of MnmE and MnmG in the presence of GDP-AIFx (cyan) is shown and 
compared with the Rg of MnmE (red) and MnmG (yellow) upon mixing with buffer. The Rg values of the a2p2 and a4p2 complexes derived from static 
SAXS measurements are indicated at the side of the Y-axis. These data show that the formation of the a4p2 complex occurs very fast, within the dead time 
of the measurement (<20 ms). Error bars show the standard deviation of three measurements. 



Our SAXS analysis corroborates this and clearly distin- 
guishes between these two conformations (Figure 1). We 
moreover observe that this transition does not only include 
the closure of the G domains but is also coupled to a move- 
ment of the helical domains toward the G domains. Our 
SAXS data further clearly support that in a first step of the 
tRNA modification cycle (Figure 7, step 1) MnmE binds to 
MnmG in an asymmetric manner, leaving one subunit va- 
cant on each MnmE and MnmG dimer (Figure 3). This ob- 
servation suggests a form of negative cooperativity between 
the subunits of MnmE and MnmG, the nature of which is at 
this point not obvious from the available crystal structures. 
However, our model for the MnmEG a2(32 complex is fur- 
ther supported by mutagenesis studies and by the disrup- 
tion of the complex by the binding of Nb_MnmG_l , which 
shares its MnmG interaction interface with MnmE. More- 
over, the same binding asymmetry is observed in the bind- 



ing of tRNA to MnmG, where only one tRNA molecule 
seems to bind to the MnmG homodimer (Figure 2). When 
GTP binds to MnmE in the complex, the conforma- 
tional changes in MnmE seem to induce allosteric changes 
on MnmG, since it promotes the binding of a second MnmE 
dimer on the opposite side of MnmG, resulting in an a4(32 
complex even at low protein concentrations (Figure 7, step 
3). The formation of this larger complex in the GTP-bound 
state is supported by high-resolution gel filtration experi- 
ments, SEC-MALS and SAXS measurements (Figure 5). In 
the next step, each MnmG monomer in the a4p2 complex 
can bind one tRNA leading to an aA^lyl complex (Fig- 
ure 7, step 4 and Supplementary Figure SIO). During (or 
prior to) GTP hydrolysis the tRNA will be modified at the 
wobble position and concomitant with GTP hydrolysis the 
large complex will dissociate again to an a2(32 form (Figure 
7, step 5). This link between GTP hydrolysis and the con- 
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Figure 7. Model of MnmEG complex formation during the tRNA modification cycle (MnmE is colored red and pink; MnmG blue and purple). (1) 
MnmE is binding to MnmG in a non-symmetric manner witli initial closure of the G domains. (2) Upon GTP binding the G domains of MnmE undergo further 
closing, with a concomitant movement of the helical domain, leading to a structural change in MnmG, (3) promoting the binding of a second MnmE dimer on 
the other MnmG subunit. (4) In the following step two tRNA molecules can bind mainly via interactions with MnmG and (5) upon GTP hydrolysis the tRNA 
becomes modified (indicated in green), the complex dissembles to its a2|i2 form and the tRNA leaves the complex. 



version from the a4(32 to the a2p2 complex is clearly shown 
in Figure 6. Considering the relatively low affinity between 
MnmE and MnmG {Ku 3-7 \iM in the presence of FAD) 
we cannot exclude that the a2p2 complex dissociates fur- 
ther into the MnmE and MnmG homodimers. Conversely, 
considering that (i) the cellular GTP concentration [1 mM- 
1.6 mM; (50,51)] is about 10 times higher than the GDP 
concentration, (ii) the GTP binding and MnmEG complex 
formation is fast and (iii) the dissociation of the G domains 
is rather slow (26), we can presume that MnmE will be for 
a considerable fraction in the 'closed' state in the cell and 
that hence the MnmEG complex will reside for a significant 
fraction in the a4p2 form. 

To unravel the physiological relevance and advantages of 
these complexities in the MnmEG-catalyzed tRNA mod- 
ification cycle, further detailed studies as well as high- 
resolution crystal structures of these protein and protein- 
tRNA complexes are required. This will also help to fur- 
ther understand the detailed molecular mechanisms behind 
the GTP-induced conformational changes promoting the 
oligomerization. 

It should be noted in this regard that MnmE is not the 
only G protein that undergoes an oligomerization step dur- 
ing its functional cycle. Dynamin and DNMIL (dynamin 
1-like protein) are, like MnmE, G proteins belonging to 
the class of G proteins activated by nucleotide-dependent 



dimerization (GADs) (20). Both proteins oligomerize via 
their a-helical stalks to form a large helical arrangement 
surrounding the 'neck' of a membrane vesicle, and dimer- 
ization of the G domains subsequently links neighboring 
turns of the dynamin helix (52-55). In contrast to these lat- 
ter two G proteins, our data suggest that MnmE does not 
oligomerize in itself via its a-helical bundle domains. Con- 
versely, two MnmE molecules use their a-helical bundles to 
bind one MnmG molecule as a connector. In the case of dy- 
namin, the oligomerization is necessary for its function in 
membrane scission during endocytosis (52). The exact role 
of the oligomerization of the MnmEG complex during the 
tRNA modification cycle remains to be explored. 
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Supplementary Data are available at NAR Online. 
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